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Summary 


A case study of a spacecraft having flexible solar arrays is presented. A 
stationkeeping attitude control mode employing both earth and rate gyro reference signals 
and a flexible vehicle dynamics modelling and implementation is discussed. 

The control system is designed to achieve both pointing accuracy and structural 
mode stability during stationkeeping maneuvers. Reduction of structural mode interactions 
over the entire mode duration is presented. The control mode employing a <Lscrete-time 
observer structure is described to show the convergence of the spacecraft attitude transients 

during AV thrusting maneuvers without pre-loading thrusting bias to the on-board control 
processor. The simulation performance using the three-axis, body-stabilized nonlinear 
dynamics is provided. 

The details of a five-body nonlinear dynamics model are discussed. The spacecraft 
is modelled as a central rigid body having cantilevered flexible antennas, a parr of flexible 
articulated solar arrays, and two gimballed momentum wheels. The vehicle is free to 
undergo unrestricted rotations and translations relative to inertial space. A direct 
implementation of the equations of motion will be compared to an indirect implementation 
that uses a symbolic manipulation software to generate rigid body equations, 
generalization of this approach to this class of flexible vehicles will be provided. 


1. Introduction 

Three-axis body stabilized spacecrafts having solar wings with significant structural 
flexibility may exhibit rigid-flex coupling effects during a typical stationkeeping maneuver. 
One of the primary concerns for the design of three-axis stabilized spacecraft is the 
structural mode interaction with the attitude control system. In addition the dynamic 
analysis and the control performance evaluation are sensitive to the ngid-tlex modelling 
accuracy. This paper presents a case study on the design, analysis and digital simulation of 
a microprocessor-based stationkeeping control system of a 3-panel communication 
spacecraft using thrusters as control actuators. It discusses the control system to achieve 
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modal stabilization, the dynamic model development and validation, and the technique for 
closed-loop digital simulation. 

The vehicle under discussion is a three-axis stabilized geosynchronous 
communication satellite. The vehicle is powered by photovoltaic solar arrays and are 
controlled in its on-orbit operation by a combination of momentum wheels, reaction control 
thrusters, and solar array motion. This satellite is modelled as a central rigid body having 
cantilevered flexible antennas, a pair of flexible articulated solar arrays, and two gimballed 
momentum wheels. The vehicle is free to undergo unrestricted rotations and translations 
relative to inertial space. The solar arrays rotate relative to the central body in response to 
the action of control torques. The momentum wheels are assumed to be controlled such 
that their motions relative to the central body are prescribed In the deployed configuration, 
the solar arrays contain 75% of the total inertia of the satellite with only 7% of the total 
mass. Each solar array has 4 cantilevered frequencies below 1 Hz. 

The solar wing flexibility is fully coupled into the body roll and yaw dynamics 
because the flexible solar wings are fixed about the roll and yaw axes of the central body, 
while articulating about its pitch axis. The pitch coupling depends upon the nonlinearity of 
the solar wing drive and its friction characteristics. The control bandwidth of the on-orbit 
normal mode is usually designed at a frequency well below the first structural mode so that 
the solar wing flexibility does not interact seriously with the normal mode controller. 
However, a relative high control bandwidth is needed to maintain pointing accuracy in the 
presence of a large thrusting disturbance. The disturbance torques are primarily induced by 
the offset of spacecraft center of gravity (CG) from the geometric or pressure center of 
maneuvering jets as well as the thrust mismatch. The sensed spacecraft flexible dynamics 
interact with the stationkeeping controllers, which may result in structural mode instability 
at high loop gains. 

.During the thrusting maneuver, the dominant modes coupling in the stationkeeping 
control are phase stabilized using the lead inherent in the sensed gyro rates together with the 
phase-lead notch filters, while the non-dominant modes at higher frequencies are gain- 
stabilized. After the thrusting maneuver, any residual rates must be nulled by an order of 
magnitude in preparation for a smooth transfer back to the normal control mode. This 
paper will focus on the stationkeeping control during the thrusting maneuver only. The 
thrusting bias about each control axis is estimated such that a fast convergence on the bias 
estimates can be achieved without an open-loop torque pre-bias by ground command, 
although this feature is also included in the design. The control loop bandwidths are 
designed as high as possible to meet pointing requirements, while still achieving adequate 
modal stability. Simulation results demonstrate stationkeeping control performance from a 
typical 2-jet (5 lbf each) south maneuver under the worst case maneuvering conditions are 
presented. 

To verify the accuracy of the dynamics model, two approaches are taken. In the 
first approach, the direct approach, the equations of motion for the vehicle are derived from 
fundamental momentum principles. The flexible appendages are modelled with 
conventional lumped mass model employing stiffness matrices to characterize the internal 
energy. Transformations to appendage modal coordinates are made and a reduction in the 
number of elastic degrees of freedom is achieved through their truncation. Simplifying 
assumptions are made regarding the magnitude of certain nonlinear kinematic terms based 
on operational considerations. The final set of governing equations are coded in a first 
order form suitable for numerical integration. 

In the second approach, the indirect approach, an unconventional method is 
employed. The rigid portion of the equations and the code is obtained from a symbolic 
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numerically compared through a set of chosen open loop comparison tests. 


2. Performance Requirements and Control System Descriptions 


Performance Requirements 


The spacecraft under study retires a, leas. 400 

”"SstT«c°sho“ bS Sch" Factors Stng pointing accuracy during the 
stationkeeping maneuver are earih sensor noise, rate gyro .noise gyro ra.el >tas «ttmat.on 

pro^itemoriom'nVxTb^ cf^T^win^^TrWector^ thnlst mismatch, thruster 

oSiTatTonSn on “ac^lfl aS.ude and provide stability range in the presence of 
structural mode frequency uncertainty. 


Control Algorithms 

The stationkeeping control is executed through special contto ^’l°p t ^y ha t t 
when the spacecraft is in Stationkeeping Mode. 

variation in thrust pulse amplitude, various thrusting disturi^ance torques a^uttheconu^ 
axes mav be induced when the maneuver thrusters are activated. Stationkeeping Moae 
provides thruster control for a maneuver execution and au tono ^ u ^Jl ltU ^ n C °J e contro i in 
body transient errors and maintain pointing accuracy during man . ^ 

roll and pitch axes is achieved using earth referenced signals and rate in t^ 
with 5 lbf thrusters. A rate integrating gyro is used as rate references in ya 

Figure 2.1 shows the functional block diagram of thej-axis ^tjonkee^ng auitude 

comm^nd^enera^ionhusi^g e^h^^et^c^p^h and Ute^rawTate 

o Vr0 drift can be minimized by an on-boaid software calibration to the raw rate 

measurements. The yaw control logic in Stationkeep.ng '^"“^‘“^ba^id on 
nitrh channels excem for the yaw gyro bias estimate, which is constant oaseu u 

premaneuver calibration estimate. Thf technique for rate gyro caltbrauon w,U no. be 
discussed in this paper. 
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The earth referenced pitch and roll position errors are sensed at 2.0345 Hz, i.e. 
every 30 real time interrupts (RTIs) at 16.384 msec per RTI, to supply position references 
for the roll and pitch gyro calibration. The angle estimate is obtained by adding the gyro 
bias estimate to the raw rate measurement, and integrating the resulting rate. For the yaw 
control loop, it does not perform the yaw gyro calibration during maneuver, since no yaw 
position position references are available in Stationkeeping Mode. Instead, yaw gyro bias 
is estimated in the normal mode prior to maneuver, and its value is held throughout the 
maneuver. 

The gyro referenced rate measurements with calibration are sensed every 8 RTIs to 
supply rate references to a third-order attitude estimator for each axis. Each attitude 
estimator performs two functions. First, it integrates the gyro rate (after correction for 
bias) to obtain a position estimate. Both the roll and pitch gyro calibration along with 
integrations operate during the premaneuver gyro calibration period as well as throughout 
the maneuver. Second, it estimates the spacecraft angular rate and acceleration bias about 
the respective axis. The roll and yaw attitude estimators also include the effect of roll-yaw 
coupling due to spinning wheel momentum. 

A proportional controller is employed for each axis to determine the control 
acceleration commands based upon the position, rate and acceleration bias estimates. The 
control acceleration commands are held constant over each control sample period. The on- 
board optimal thruster selection (OTS) logic selects available thrusters and determines 
necessary thruster on/off command duration to valve drivers based upon the minimum fuel 
consumption. The selected 5 lbf thrusters are turned on/off for commanded durations to 
deliver the control momentum equivalent to the commands, and achieve attitude corrections 
during maneuvers. The detailed technique for conversion of the control acceleration 
commands to thruster commands will not be discussed further. Although the control 
sample period is designed to be commandable, it must be selected to meet the needs of the 
control processor thruput and avoid structural mode instabilities as well. 


3. Flexible Spacecraft Model Descriptions 

The mechanical idealization of the satellite is illustrated in Figure 3.1. The model 
consists of a central body, which is considered to be rigid and to which are mounted a set 
of reaction control thrusters. Cantilevered to the central body are a complement of 
structurally flexible antennas. Two distinct, structurally flexible solar arrays are hinge 
connected to the central body. The arrays can rotate independently about parallel drive axes 
in response to control torques, which are assumed to be known functions of time. Two 
independent, identical, rigid axisymmetric, variable-speed momentum wheels are mounted 
to the central body through two-axis gimbal mechanisms. The complete motions of the 
wheels relative to the central body are assumed to be prescribed functions of time. 

The communications antennas of the actual satellite are capable of limited 
articulation relative to the central body. However, because these rotations are small in both 
magnitude and rate, their influence on the vehicle's overall attitude dynamics was deemed 
negligible, and these degrees of freedom were not included in the model. While the solar 
arrays will be virtually identical under nominal circumstances, they are treated as 
structurally distinct to accommodate more general conditions. 
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The deformations of the antenna and solar array structures are assumed to be linear 
elastic in character and small in magnitude. The respective appendages are modelled as 
collections of point masses interconnected by massless elastic structure. Stiffness matrices 
are used to define the elastic restoring forces acting internal to these assemblies. 
Ultimately, modal coordinate transformations are introduced for each appendage and the 
final motion equations are cast in terms of truncated sets of those variables. 

Cantilevered frequencies of an individual solar array and of an assembled antennas 
are provided in Table 3.1. 

Table 3.1. Characteristic Cantilevered Frequencies of the Appendages. 


Mode # 

Solar Array 
Frequency (Hz) 

Antennas 
Frequency (Hz) 

1 

.118 

1.586 

2 

.355 

1.792 

3 

.705 

1.953 

4 

.835 

2.043 

5 

1.825 

4.235 

6 

2.725 

4.867 

7 

3.167 

5.323 

8 

4.914 

5.893 

9 

6.161 

13.366 

10 

6.905 

19.558 
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4. Control Design and Analysis 


Linearized Spacecraft Open Loop Dynamics 

Let { b } be the spacecraft body frame, (s)w be the solar wing frame and C w be the 
direction cosine matrix at the wing-to-body angle a such that (s} w = C w {b}, where w = n 
(North wing) or s (South wing). Figure 4.1 is a geometry showing the spacecraft with the 
thrusting forces and torques and the disturbance created by thruster plume impingement; in 
which we assume that point i is the pressure center on North wing where the resultant 
plume force vector applies. 



Figure 4.1 A Geometry showing External Thrusting Forces/Torques 
and Plume Disturbance 


The linearized body-stabilized spacecraft hybrid dynamics with articulated solar wings are 
given in Equations (4.1) - (4.6). 


m x + P n r| n + P s r) s = F l + R 
L/C cil + In e 2 cbn + Is e 2 ci) s + On iin + Q s rj s 
4 I n CO + 15 (On + e? Qg T\ n = Tf* 

4 I s 0) + 1*2 co s + 4 Q£ ri, = T| wd 


(S/C translational motion) (4.1) 

- T 1 + TP (S/C rotational motion) (4.2) 
(North wing pitch dynamics) (4.3) 

(South wing pitch dynamics) (4.4) 
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2 t C n F 

Pi X + Qj 0) + (e^ QS) T ton + Tin + 2CA Tin + A Tin = 

L a ly 

(North wing flex dyn.) (4.5) 

pj x + Ql co + (e£ C§) T cb s + ris + 2£A q s +Ati s = 0 (South wing flex dyn.) (4.6) 
where 

m = total S/C mass 

Is/ C = S/C mass inertias 

x = S/C translational position vector 

co = S/C angular rate vector 

I n ,I s = pitch inertias of wings about their interface points 
co n ,co s = relative pitch angular rates of wings 
rin,T)s = modal variables of wings 

P° ( Pf = rigid-flex translational coupling matrices of wings about their 
interface points 

P n P s = rigid-flex translational coupling matrices of wings about the 
S/C CG 

Q£, Cg = rigid-flex rotational coupling matrices of wings about their 
interface points 

On, Qs = rigid- flex rotational coupling matrices of wings about the 
S/C CG 

jswd -p|wd _ so)ar win g tor q Ues 

FP.ATP = plume induced force and torque about pressure center 

(J> = mode shape at point i of North wing (a nx6 matrix with 3 translational 

and 3 rotational deformation, where n is the number of modes in 
concern) 

A = cantilever mode frequencies for each wing 

£ = structural damping factor 

e 2 =[0 1 0] T 


The model above is with articulated solar wings driven by the wing torques about 
their hinge axes. To fully include the flexibility of wings, the dynamic inertias of each 
wing referenced to its interface point must contain at least 99% of the roll or yaw moment 
of inertia about the same point, or, the pitch inertia excluding yoke. The flexible reflectors 
have relatively small dynamic inertias compared to the total spacecraft mass inertias; the 
structural mode control interaction is negligible. Therefore, the reflectors are considered to 
be rigid and included as a part of rigid central body in stability analysis For the pitch 
dynamics of wings in the ’’constraint" state (see Design Considerations ^ow for further 
discussion), the solar wing angular acceleration terms in Equations (4.1), (4.2), (4.4) ana 
(4.5) may be eliminated from these equations, and the resulting model represents the 
linearized spacecraft dynamics with non-articulated wings. 
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Control Design Model of Rate Loops 

The control system has an outer loop (i.e. position loop) and an inner loop (i.e. rate 
loop). The position loop is designed at very low bandwidth with its gain crossover well 
below the structural modes and consequently has generous gain margin (> 40 dB) on all 
flexible modes. As for the modal stability, the rate loops are the primary concerns. The 
design model of rate control loops in Stationkeeping Mode is given below. 


Rate Gyro Dynamics 


ef 


a, 


= 1 

s s + a 


g 


g 


to i 


(i= 1,2,3) 


Gyro Processors co®(n+l) - [9f(n+l) - 0f(n)J / AT 

~g ~g 

Attitude Estimators 0i( n+ l) =0i(n) + AT (D?(n+1) 

cof(n+l) =cof(n) + Aco'(n) + ATd;(n) + Aco^(n) 
cof(n+l) = tof(n+l) + K r [ojf(n+l) - cof(n+l)] 


dj(n+l) = dj(n) + K<j [cof(n+l) - cof(n+l)] 

Aco^(n) = -(HjAT/Ii) ti^(n) (i=l) 

= 0 (i=2) 

= (H t AT/I 3 ) co?(n) (i=3) 

Aco|(n) =ATaP f (n-l) 

Proportional Controllers = *Np(z) (c p 9i + C r O)f) - d,(z) 


N p (z) 


n2Z 2 + njz + np 
z 2 + d]Z + do 


(Phase Lead Notch Filter) 


Control Transport Delay 


a T = ^otf f (n-p) if n AT + x d < t < n AT + t d + x c , p = (x d + .5 x c )/AT 

ic 

— d otherwise 

where 

0j = spacecraft angular position about body axis i [deg] 

0)i = spacecraft angular rate about body axis i [deg/sec] 

®i = estimate of 0j [deg] 
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oaf = rate gyro measurement with calibration [deg/sec] 
cof = estimate of coj [deg/sec] 

di = acceleration bias estimate about body axis i [deg/sec 2 ] 

ef = position signal output from the rate integration gyro [deg] 

RT1 = real time interrupt = 16.384 msec 

of = filtered control acceleration commands 

aj = actual control acceleration acting about body axis i over x c [deg/sec 2 ] 

Aco] = feedforward rate changes [deg/sec] 

AT = control sample period [sec] 

Td = control pulsing delay [sec] 

Xc = actual control pulsewidth [sec] 
n = control sampling time in Stationkeeping Mode 
K r = rate estimation gain [(deg/sec)/(deg/sec)] 

K<j = acceleration bias estimation gain [(deg/sec 2 )/(deg/sec)] 

Cp = spacecraft position control gain in Stationkeeping Mode [sec ] 

C r = spacecraft rate control gain in Stationkeeping Mode [sec -1 ] 

N p (z) = phase-lead notch filter, 
ag = gyro servo bandwidth [rad/sec] 

[ I = spacecraft moment of inertias about body axis i [slug-ft ] 

Hj = total spacecraft angular momentum about pitch axis 


It employs a discrete predictor-corrector algorithm to estimate the spacecraft rate and 
acceleration bias and the integration of the rate measurement(after correction f or^robias) 
for position estimate. A discrete phase-lead notch filter is added in senes with the attitude 
control acceleration command (i.e. the proportional controller excluding the acceleration 
bias control term, which is added to the filtered attitude control acceleration conunand) to 
provide additional phase lead for modal stabilization. The spacecraft dynamics used in the 
design are represented by the linearized hybrid dynamic model. The rate gyro modefi 
which has a first order servo of 8 Hz bandwidth, output positional signal, and the gyro 
processor determines the rate based upon the position change over OTC * 

period. The equivalent transfer function of the design model will not be given in this 
paper. We will discuss various design concerns which are related to transient performance 
and structural mode stability in general, and the design philosophy in achieving the goals. 


Sensitivity To Modal Parameters 

The stability of the structural modes selected for baseline design may be affected by 
both the structural frequency uncertainty and the structural damping. A structural damping 
ratio of 0.0025-0.005 is added to the hybrid dynamic model. The P and Q matrices defined 
earlier are essentially the diagonal matrix elements of the translanonal and rotational ngi - 
flex coupling matrix B given in Equation (6.9). The coupling matrix selected for the 
baseline design is derived based upon the the spacecraft on orbit nominal configuration 
with fully deployed wings whose z-axis is directed to the Earth. By knowing thelocation 
of the interface point relative to the spacecraft CG and the wing orientation, it can be proved 
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that conversion of the coupling parameters to about the spacecraft CG is accomplished 
through the relations: 


p n = Cj PS , On = Cjf OS + Rn Pn 

P s = Cl P° , Qs = Cl Q° + R s P s 


(4.7) 

(4.8) 


where ~ denotes the skew symmetric matrix operator which achieves a vector cross 
product. P and Q will change when the spacecraft is no longer in nominal configuration 
under the following two conditions. First, P and Q vary as the wings rotate about their 
hinge axes; and second, due to the nonlinear characteristics of panel hinge stiffness. The 
solar wings could be at any orientation with 16 deg or less wing separation angle Rotation 
of wings will primarily affect P and Q about the spacecraft body roll and pitch axes, and 
almost no change about the body pitch axis in the presence of the symmetric north/south 
wings. Furthermore, when the east or west thrusters fire during an east/west maneuver 
with non-zero wing angles, the panel hinge loads as induced primarily by the linear 
acceleration of the spacecraft along the the panel z-axis as well as the flexibility of wings 
may exceed the spring preload such that the panel stiffness will drop from its hardstop 
region, where the nominal P and Q are derived, to the deadband region, where a soft panel 
stiffness is present. The worst case panel hinge loads result when the wing is at 90 deg 
orientation during an east/west maneuver. When this occurs, P° and Q° about the interface 
point of the wing will vary about all three axes. Both the roll and yaw control axes must be 
designed to stabilize all possible structural mode frequencies which may result from the 
rotation of wings and the nonlinear characteristics of panel stiffness. 


Spacecraft CG Uncertainty 

4Tie CG offset of the spacecraft from the pressure center of maneuver thrusters will 
result in a thrusting disturbance about the control axis. If the actual CG offset was 
predicted to a 100% accuracy by the OTS in advance, the thrusting disturbance would be 
self-compensated with the selected thrusters - the unique feature of the OTS. The CG 
uncertainty of the spacecraft has a major impact on the maneuver transients. The concerns 
are in two areas: (i) the CG uncertainty of the spacecraft while on station at a steady state 
condition, which is primarily caused by the tank misalignment and the possible imbalance 
of the dry spacecraft and (ii) the CG uncertainty due to propellant motion during a 
maneuver. The estimation errors on the acceleration bias may result in an initial transient 
about each axis that exceeds steady state pointing. The acceleration bias estimation gain 
(Kd) must be designed to minimize the maneuver initialization transient, and to avoid 
excitation of structural modes, while still being able to track disturbance. The transient 
errors can also be improved with an initialization of the acceleration bias estimates to the 
steady state values recorded from the last maneuver. 


Firing Thrusters on Flexure 

The wing flexibility may be fully excited in the steady state condition from a long 
maneuver bum. Also, when switching the control logic to further null the rigid body 
residual rates at the completion of a maneuver, an instant loss of thrusting forces could 
yield a significant response to the already excited flexible dynamics with the magnitude 
exceeding the impulse control deadband limits. To avoid firing on flexure the rate 
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estimation gain (K r ) in Stationkeeping Mode must be as low as possible to filter the sensed 
spacecraft flexible dynamics, so that the proportional controllers determine the required 
control momentum based upon the estimated rigid body dynamics. The rate gain in the 
range 0 < K r < 0.5 (deg/sec)/(deg/sec) meet the requirement, while still being able to track 
the spacecraft dynamics to a degree of accuracy by feeding forward the commanded rate 
changes to the estimators from the OTS. 


Control Transport Delay 

Due to onboard control software processor speed, a computational delay on rate 
change command processing in OTS is induced. Such a delay together with the phase 
delay induced by the gyro rate signal processing and the control pulsewidth delivery induce 
a phase lag to each control loop, which affects the stability of structural modes. To 
compensate for the loss of phase due to the control transport delay, a phase-lead notch filter 
is employed in Stationkeeping Mode to provide each control loop with an additional phase 
lead The design philosophy is to set the modulation frequency so that the half sample 
frequency is well above the dominant modes to ensure phase stabilization of these modes 
with additional phase lead produced by the notch filter. Because the notch fiber is not wide 
enough to provide all modes with sufficient phase lead, the modes near the half sample rate 
are gain-stabilized. Phase stabilization of the dominant modes simply means that the 
control loop will generate a stabilizing feedback signal to that mode. This is to be 
contrasted to gain stabilization wherein the non-dominant modes rely upon the structural 
damping of the spacecraft to provide enough damping to overcome any slight destabilizing 

effects. 


Effect of SWD Deadband on Torsional Mode Stability 

Due to the solar wing drive (SWD) backlash, the flexible pitch dynamics of wings 
about their interface points may appear in one of the three states: free-free , constraint 

and one between these two states, depending upon whether the wings are inside or outside 
the deadband and the magnitude of friction. The SWD has a 0.5° deadband. When the 
wings are inside the deadband and the magnitude of friction is insufficient to overcome the 
wing relative motion, the central body pitch dynamics are then disturbed by the load 
torques with a phase-shifted bang-bang profile whose magnitude is equal to Coulomb 
friction It is very complex to analyze the torsional mode stability with such a profile, lhe 
best way to examine the performance is through the simulation by actually including the 
SWD When the wings are inside the deadband but the friction is negligible so that the 
flexible dynamics of each wing are "free-free" about its pitch axis, the central body pitch 
dynamics are then fully decoupled from the wing relative motion. The central body under 
this state is considered to be rigid. The stabilities of the free-free modes will not become 
a problem as long as the free-free motion remains inside the deadband. Still, when inside 
the deadband but the relative motion of the wing is locked up by friction, or, the wing-to- 
body rates are large enough to break the friction and the gear teeth are recontacted to the 
SWD shaft, the pitch flexible dynamics of wing are now considered to be constrained 
about its pitch axis; that is, the SWD will output load torques in absence of stepping 
commands, to drive the wings to prevent them from moving about the pitch axis of the 
central body. From this viewpoint, the wings are also fixed to the central body about its 
pitch axis and the flexible dynamics of wings are fully coupled into the body pitch axis 
through the load torques. 
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Plume Impingement Effects 


Forces and torques created by thruster plume impingement on spacecraft 
appendages (north wing and east/west antennas) may result in a net change in the overall 
contto! torques and an unmodelled excitation of the structural modes. From the spacecraft 
stability viewpoint, the principle concern is that the plume disturbance shall not add phases 
to erode the phase margins of the phase stabilization modes nor magnify amplitudes to 
degrade the gain margins of the gain stabilization modes. 

Since r is the pressure center of plume forces, the plume induced torques about the 
deformed S/C CG is 


T p = (R n + r ) x F p + AT P 

AT p = T p - (R n + r ) x F p = {s} t {c„TP - (QR n + r)(C n FP)) 
ATP = C n T p - (QR n + ?tc n FP) 


Referring to Equation (4.5), the modal excitation of the north wing is induced by (i) the 
spacecraft and the north wing motion driven by the terms 

e n = Pi x + Qj © + (ej C£ ) T o>n 
and (ii) the plume disturbance 


e p =<I> T [(C n FP) T (ATP) T ] T 

acting about point i, the assumed pressure center of plume forces. 

In the presence of the plume disturbance eP acting on the North wing, the open 

loop dynamics transfer function g^s/T 1 (s) from an impulse response can be derived from 
the single axis hybrid dynamic model with one mode only as 

to(s) _ | l-pi s 2 + 2 ^]Xis + X] 

TW ' >*/c 1-1*2 S 2 + 2 W2S + x § * 

where 

Pi=ePq n /T‘ 

M^2 = 2 / I s / C 

Cl =c/^rr 
C 2 = C / V 1 -H2 
h = X / V 1-M-, 

h =x/vr^ 
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The parameters eP and q n are the plume disturbance and the rigid-flex rotational coupling 

term associated with the structural mode of frequency X and modal damping T 1 is the 

impulse thrusting torque about the single axis in concern. If Hi * 0 (i.e. eP * 0), the plume 
disturbance will perturb the zeros of the transfer function above from their nominal 

locations. If the perturbed zeros move toward the poles (i.e. Hi > 0) , then it improves the 
margin of the phase- stabilized mode, producing less rigid-flex coupling. On the contrary, 

if Hi < 0, the zeros move away from the poles, and the loop gain is magnified by a factor 
of (1- Hi )'■ both of these factors will erode the margin of any gain-stabilized mode. 

The geometry of the thrusters and solar wings are such that, in fundamental modes 
(those with no inflection points), the plume impingement coupling and the ngid-flex 
dynamic coupling act in phase with one another to excite a mode. In other words, plume 
impingement acts to amplify modal excitation already present due to ngid-flex dynamic 

coupling. This implies that Hi < 0 for fundamental modes. In this case, the zeros of the 
transfer function move even further from the poles, exacerbating the flexible dynamics 

coupling problem for these modes. Thus if Hi < 0> it is desirable to have the magnitude of 

Hi as small as possible: Hi 1 « 1 H2 I is goodness. In this case study, Hi = -0.0388 for the 
first out-of-plane mode at 0.1 185 Hz. Fortunately, values are small in comparison to the 

corresponding H2 = 0.7037, indicating that plume impingement is not a dominant effect. 
One measure is the zero/pole frequency ratio: with no plume impingement 


= YTH 2 = Y1-.7037 = .544 
X2 


whereas with plume impingement 


X) _ 
X2 



1-.7037 

1+.0388 


.534 


This is a change of only 1 .87% which is small relative to the 5% or 10% accuracy to which 
H2 is known to begin with (A 1.58% change in the value of H2 would result in the same 
change in zero/pole ratio). 

Higher frequency modes whose mode shapes include an odd number of inflection 
points between the attach point and the "point of application" of the plume impingement 

force can exhibit a positive value for Hi- 1° this case, the plume impingement force acts 
opposite to the direct rigid-flex dynamic coupling and tends to reduce modal excitation. In 
the transfer function this is reflected by the fact that the zeros move closer to the poles, thus 

tending to cancel. Should Hi ever get as big as H2> the zero would exactly cancel. An even 

higher value of Hi would reverse the phase of the modal coupling. This situation is of less 
interest, and is probably not possible with plume impingement as the excitation source. 
Since plume impingement is, in reality, a distributed force rather than a point force as 
modeled herein, its viability in exciting a higher frequency mode diminishes rapidly as the 
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number of inflection points increases. The assumption made herein that plume 
impingement force is applied at a single point loses its validity for higher frequency modes, 
therefore results should not be taken too literally for such modes. The "constrained" state 
is similar to those for roll or yaw loop. 


5. Control Loop Stability 

With the panel support cantilevered at its base about the transverse axes, but free in 
torsion, Table 5.1 characterizes the flexibility of 0 deg, 3-panel single wing in terms of the 
modal frequencies with associated dynamic inertias about the interface point of wing. The 
first twelve modes as listed contain > 99% of the total inertias of wing about each axis, 
which are sufficient to describe the flexible characteristics of wings. When the wings are 
attached to their base, the flexible dynamics appeared to the angular motion of the 
spacecraft through rigid-flex coupling have frequencies higher than that of the cantilever 
modes. The increased modal frequencies, assuming perfectly symmetric wings, are 
defined as the system modes in Table 5.1. The frequencies of the system modes will shift 
as the wings rotate about the hinge axes; 90 deg wings yield out-of-plane cantilever modes 
in the yaw axis and in-plane cantilever modes in the roll axis. For the transfer function of 
the coiTesponding open loop dynamics, the dominant modes have wider pole/zero 
separations. One of the design goals is to stabilize the system modes under any wing 
orientation. 


Table 5.1. Solar Wing Structural Modes 

at Zeno Degree Wing Angle 


Frequency, Hz 

Dynamic Inertia 




System 

(about interface 


Mode 

Cantilever 

Constraint 

Free-Free 

point), Kg-m 2 

Definition 

1 

0.1185 

0.2117 


44.498 

oul-of-plane 

2 

0.3547 

0.6323 


45.970 

in-plane 

3 

0.7051 

0.7631 


8.722 

out-of-plane 

4 

0.8508 

0.8582 

1.3645 

4.498 

torsional 

5 

1.8254 

1.8627 


3.550 

out-of-plane 

6 

2.8058 

2.8091 

3.1938 

1.588 

torsional 

7 

3.1668 

3.2573 


3.412 

out-of-plane 

8 

5.1320 

5.1338 

5.3698 

0.854 

torsional 

9 

6.1608 

6.1704 


0.592 

out-of-plane 

10 

8.2235 

8.2333 


0.432 

out-of-plane 

11 

8.9668 

8.9898 

9.1631 

0.643 

torsional 

12 

9.4732 

9.4882 

9.7306 

0.524 

torsional 
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Figure 5.1 shows the discrete-time Bode plots and Nichols chart of the 8 R 
Stationkeeping Mode spacecraft roll rate control loop with no structural filter or conffol 
transport delay. The control bandwidth was designed to limit transient errors to within 0.1 
deg in the presence of a 1.5 inches spacecraft CG offset along the z-axis. The first five 
out-of-plane modes at nominal frequencies, 0.5% structural damping and 0 deg wing ang 
were included. Using gyro references, both the first (mode 1) and second out-of-plane 
mode (mode 3) are phase-stabilized with about 70 deg and 30 deg phase margi , 
respectively, and the remaining out-of-plane modes (5.7, 9 & 10) are gain-stabilized with 
at least 22 dB gain margin. The control design provides a 6.5.1 ratio to the separatio 
between the zero gain crossing frequency and the pole of the first structural Figure 

5.2 shows the same design without a structural filter, but with a 2 RT1 control ^^port 
delay. The transport delay effect can be seen on the structural modes greater than 0.35 H , 
to which the phase lag induced by transport delay was added, yielding almost no P |ase 
margin on the second out-of-plane mode. The linear design was then improved with a 
phase-lead notch filter, which has a unit gain in the low frequency range and a maximum of 
68 deg phase lead at the notch frequency of 1.209 Hz. Figure 5 3 shows that the phase 
margin of the second out-of-plane mode was increased up to 36 deg with the phase-lead 
notch filter. The notch frequency was carefully selected to ensure that all phase-stabilized 
modes will remain in the phase stabilization region in the presence of 100% frequency 
increase as shown in Figure 5.4. Although a 100% frequency increase is allowed before 
mode 3 loses its phase stabilization characteristics, this mode is also gain stabilized once its 
frequency increases from the nominal. Also, the third out-of-plane mode (mode 5) is gam 
stabilized with 13 dB margin at its nominal frequency (1.8627 Hz system mode). 
Decreasing the frequency of this mode immediately leads it to the phase stabilization region, 
while the gain stabilization characteristics are still retained. This mode reaches adequate 
phase margin before the gain stabilization characteristics vanish at more than 50 h 
frequency drop as shown in Figure 5.5. 

The pitch transient during south maneuvers is affected by thrusting disturbance, 
primarily induced by both the thruster cant angles and the thrust mismatch. The net pitch 
disturbance is estimated to be 0.4 ft-lb, which requires the 8 msec minimum control 
pulsewidth to be fired at a rate of 1.667 Hz at which the half sample rate is nearly equal to 
the first torsional mode. The pitch loop has rigid response m south stationkeepingbecause 
the effect of the torsional mode is insignificant about the half sample wntrol rate. The pitch 
transient during east/west maneuvers is primarily affected by the S/C CG yaw offset from 
the pressure center of the maneuver thrusters. With a 5.9 inches yaw offset, it requires a 
1 5 RT1 control pulsewidth to be fired every modulation period. The linear frequency 
analysis of the spacecraft pitch rate loop in Stationkeeping Mode when the pitch wing 
dynamics with SWD are in the "constraint" state is shown in Figure 5.6. The design with 
the same phase-lead notch filter as applied to the roll and yaw loops was based on 2 RTIs 
control pulsewidth, 2 RTIs control transport delay and 8 RTIs modulation period. The first 
torsional mode (mode 4) is both phase and gain stabilized at its nominal frequency with 3>55 
sup o$ phase margin and 20 dB gain margin. Either dropping or increasing the frequency 
produces no impact on stability at all. 

Figure 5.7 shows the stability of the first in-plane mode (mode 2) at nominal 
frequency and with the phase-lead notch filter. It is also phase-stabilized with 46 deg phase 
margin. This mode remains in the phase stabilization region even with a 100% frequency 
increase while still having 23° phase margin as shown in Figure 5.8. 
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Figure 5. 1 Discrete Bode Plots and Nichols Chart of 
Roll Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, Nominal Frequencies : without Structural 
Filter and Zero Control Transport Delay 


Ruse, deg 


Figure 5.2 Discrete Bode Plots and Nichols Chart of 
Roll Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, Nominal Frequencies : without Structural 
Filter, but with 2 RTIs Control Transport Delay 
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Figure 5.3 Discrete Bode Plots and Nichols Chart of 
Roll Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, Nominal Frequencies ' with Phase-Lead 
Notch Filter and 2 RTIs Control Transport Delay 




Figure 5.4 Discrete Bode Plots and Nichols Chart of 
Roll Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping. 100% Frequency Increase : with Phase-Lead 

Notch Filter and 2 RTIs Control Transport Delay 
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Figure 5.5 Discrete Bode Plots and Nichols Chart of 
Roll Rale Loop in Stationkeeping Mode at 0.5% Structural 
Damping, 50% Frequency drop : with Phase-Lead 
Notch Filter and 2 RTIs Control Transport Delay 




Figure 5.6 Discrete Bode Plots and Nichols Chart of 
Pitch Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, Nominal Frequencies : with Phase-Lead 
Notch Filter and 2 RTIs Control T ransport Delay 
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Figure 5.7 Discrete Bode Plots and Nichols Chart of 
Yaw Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, Nominal Frequencies : with Phase-Lead 

Notch Filter and 2 RTIs Control Transport Delay 



Figure 5.8 Discrete Bode Plots and Nichols Chart of 
Yaw Rate Loop in Stationkeeping Mode at 0.5% Structural 
Damping, 100% Frequency Increase : with Phase-Lead 

Notch Filter and 2 RTIs Control Transport Delay 
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6. Dynamic Model Validation and Digital Implementation 

The main objective is to consider issues dealing with flexibility in multibody 
dynamics. Multibody dynamics is differentiated from structural dynamics by its capability 
to undergo arbitral^ rigid body motion. Analysis of flexible structures are well established 
using finite element method within the context of structural dynamics. In multibody 
dynamics context, need for flexibility modelling arose in recent years as exemplified by 
large space structures and 3-axes stabilized satellites. 

To systematically address the addition of flexible domain, the virtual work principle 
is chosen as the basis for derivation. The motivation for this choice is based on intended 
discretization using the finite element method. By choosing the same basis for multibody 
dynamics and for the finite element method, extensions into nonlinear flexibility is natural 
and consistent. Other choices are readily available in the literature [1-3]. 


Virtual Work Principle 


An integral representation of the governing equations of motion of solids are 
imbedded in the virtual work principle. By deriving the multibody dynamics equations via 
the virtual work principle, a consistent treatment of flexible domain can be made. In 
practice, the flexible domain is discretized using the finite element method. The virtual 
work principle is the basis for the finite element method. The technology developed in the 
finite element method can be integrated into the flexible multibody dynamics efforts. 

The virtual work principle states 


where 



I. 


= 8W im = 8e : odV 


R = material particle position vector wrt inertia] frame 
f = force/unit volume 

p = density 

V = reference configuration 
e = strain 


o = stress 


( 6 . 1 ) 


The main advantages offered by applying the virtual work principle are twofold. 
First, the integral representation together with the virtual displacements allow domain 
decomposition between the rigid and the flexible portions of a vehicle. Second, a 
consistent formulation of a flexible multibody vehicle can be derived and assessed. 
Consistency refers to final discretization using the finite element method. Once such 
consistent derivation is made, extensions to nonlinear flexible models can be made by 
adopting techniques developed in the finite element method [14]. 
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Rigid Body with Attached Flexible Appendage 

To further explore the method presented by the virtual work principle, an idealized 
flexible spacecraft model is derived. The idealization involve representing the spacecraft as 
a rigid body with attached flexible appendage. The flexible appendage is assumed to be 
fixed to the rigid body. Articulation is not allowed. Even with this simphtying 
assumptions, a wide class of vehicles can be modelled. 

Consider the idealization shown in Figure 6.1. 



The domains, frames, and the position vectors are defined as 

Vp = flexible domain 
Vr = rigid domain 

T\ = inertial frame 
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fb = body frame 

R = inertia] frame origin to body frame origin 

- inertial frame origin to rigid body material particle 
rR = body frame origin to rigid body material particle 

s = body frame origin to rigid body center of mass 

= inertial frame origin to flexible material particle 
rF = body frame origin to reference material particle position 
T] = relative particle displacement 

c = body frame origin to vehicle center of mass 

Application of the virtual work principle to this vehicle yields 


I 5^r ■ (f - p^ R )dV + J 8 ^f ■ (f - p^jdV = f 5e : a 

'V R Jvf Jv F 


dV 


(6.2) 


with 


and 


= R +r R 

= Fb [ii + & x a + ia x rR + ix(ojxr R )] =Fb Sr 
5£r = f b [sx +50 X rR ] 

- R +r F + ti 


^F = Fb[u + £0xu + £2x(rF + H ) + 2 & x ij + & x (& x ( rp + tj )) + tj] = Fb ap 
5^f = Fb [5x +50 x ( rp + n ) + 5a] 

where 

u = velocity of body frame wrt inertial frame 

& = angular rate of body frame wrt inertial frame 

5x = virtual displacement of the body frame 
50 = virtual rotation of the body frame 

5a = virtual relative displacement 

The components, underlined, are defined with respect to the first occurrence of the frame 
definition. For example, in above definition, the components u and ££ are defined in the 
body frame. For further discussion on this notation is clearly covered in [1], Note that the 
virtual quantities are obtained through infinitesimal variation of the current equilibriated 
state. Substituting the above quantities into Equation (6.1), following three sets of 
equations, with respect to the body frame, can be derived. 


3 \ 4 - 


F = m]i + m& xii + mmx£+m&x(&><£)+2&x p d V j + I ilpdV 


(LH*(L a 


(6.3) 


T = m£X]i + m £x(QXii) + Ii!j + £Q x I 

+ 2 ( ef + n) x (a x n) p dv + ( if + n) x n p dv 

J vf vv f 


(6.4) 


/■ 

/V F 


5n T If dv = [ | 5n T p dv 


((; 


(1 


a + ( 5n pdv 

vp / 


(I; 


+ 5n a x ( ef + 


(I: 


n ) p dvj +H 8n T ax(fijx(rF+ n)) pdvj 

(I 1 --") 


+ 2(1 5h T ( & X 11 ) p dV | + | j 8q t i p d V ] + 


with 


vf 


^Vp 


(6.5) 


F = F b + I fpdV 

/v F 


I = lb + (iF+n)xfFdv 


/■ 


where 


Fb = force applied to the rigid body at body frame origin 
Tb = torque applied to the rigid body about body frame origin 

I = instantaneous vehicle inertia matrix wrt body frame origin 


Often in practice, the integral representation in Equations (6.3) - (6.5) is skipped by 
assuming the lumped mass idealization. However, the discretization of the flexible domain 
into finite element idealization stem from these equations. The lumped mass idealization is 
an extreme case. Such inconsistent assumption with finite element method may produce 
inaccurate results for crude finite element mesh. More systematic study should be made to 
assess the consequence of such assumption. 

As closure, the lumped mass idealization will be made to produce a set of equations 
that may be compared to previous derivation [4]. The lumped mass idealization takes the 
volume integral and cast it into a sum spanning the total number of nodes in a finite element 
mesh. For an arbitrary function, this idealization can be expressed as 




where 


/( a, n, n)p dv = £/( a „ ^ )m, 


3i = finite dement nodal displacement vector 
mj = corresponding lumped mass 

Adopting this idealization. Equations (6.3) - (6.5) can be reduced to 

F = mji + maixii + m^x£ + m & x (& x £ ) + 2 £ m ; & x & +£ m iOi 

' i (6.6) 

T=mcxu+m cx(axu) + Ia+ 

+ 2 Z m > ( li + 3i) x (ffl X (ji) + X m i ( Ei + 3i) X Qi 

1 1 (6.7) 

fi = mai + m i axu + niiax(r J + a,) + mj & x (aj x ( £ + ai )) 

+ 2 ITl; X ai + m; ai + S K ijaj 

J (6.8) 


where 

m = total vehicle mass 

Kjj = assembled stiffness matrix 

By interpreting the stiffness matrix as the tangent stiffness matrix, the equations are valid 
for nonlinear flexible systems. Since modal reduction generally is not possible for 
nonlinear flexible systems, the finite element nodal degrees of freedom must be used to 
represent flexible degrees of freedom. For linear flexible system, an indepth coverage of 
an alternate derivation of Equations (6.6) - (6.8) is provided in [4], 


Extending Symbolic Rigid Body Code to include Flexibility 

In the recent years symbolic manipulation software capable of generating rigid body 
code became available. Some example of such codes are SD/FAST (Symbolic Dynamics, 
Inc.) [6], AUTOLEV (OnLine Dynamics, Inc.) [8], and AUTOSIM (Univ. of Michigan) 
[9], For rigid vehicles, these tools can dramatically reduce the time spent on deriving and 
implementing the equations of motion. 

By combining the codes generated by the symbolic manipulation software with 
reduced set of "hand" derived equations addressing the flexible domain, the capability of 
these codes can be extended to flexible vehicles. A systematic method for such an 
extension is provided for a satellite class of flexible vehicles. This method will be 
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illustrated with previously derived equations of motion for a rigid body with flexible 
appendages. 

For a rigid body with flexible appendages, the resulting equations can be partitioned 

into 


A B jii R j _ | Rr + Rrf | 

b t , J 1 ii F j 1 Rf ‘ 


where 

ur = rigid degrees of freedom 

up = modal amplitude degrees of freedom 


(6.9) 


Note that the P and Q submatrices are defined in the previous section. In Equation (6.9), a 
modal reduction has been assumed. The portion of the partitioned equation generated by a 
symbolic manipulation software is 


[A] {ur} = {Rr} 


( 6 . 10 ) 


This portion is obtained by supplying the symbolic manipulation software information on 
the current configuration. In another words, the total vehicle is assumed to be rigid. The 
requirement of current configuration entails configuration update at each integration step. 
The current configuration is the reference configuration. 

The solution process follows by forming 


{ up) = (Rf) - [b t ] {ur} 


( 6 . 11 ) 


and substituting into the rigid partition to yield 


[a -BB t ] {ur) = {Rr + Rrf} - [B] {Rf} (6.12) 


In terms of actual equations of motion, by observing the structure of Equations (6.4) - 
(6.6), the necessary additional partitions can be generated by discretization and modal 
reduction of the terms 
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I/.' 


(/. 


npdv 


(if+ n)x iipdv 


[B] {ii F > 


(6.13) 


I fpdv 


2 £2 X 1 1 a p dV 

' wvp 1 

\ I [ \ 

j 

\K 1 

if \ i 

(if+ n) x f F dv 

| | 

2 1 ( ip + n) x (« x a) P dv 

1 yv F / j 


IWvp If 


(Rrf) 


(6.14) 


The flexible partition can be generated by Equation (6.5). 

Methodology presented above produces an "exact" set of equations. Standard 
assumptions such as constant vehicle center of mass and inertia together with small relative 
flexible displacements can be made as deemed plausible to reduce computational effort. 


Rigid Body with Articulated Flexible Appendages 

With assumption that an symbolic manipulation software will be used to generate 
the rigid partition of the equation of motion, only the required matrices for the articulated 
flexible domain will be documented. The idealized articulated flexible appendage is shown 
in Figure 6.2. 
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The quantities are defined as 

Vf = flexible domain 

T\ - inertial frame 

fb = body frame 

Fk = appendage frame imbedded in the yoke body 
R = inertial frame origin to body frame origin 

A = body frame origin to appendage frame origin 

^F = inertial frame origin to flexible material particle 

rp = appendage frame origin to reference material particle position 

r| = relative particle displacement 

c = body frame origin to vehicle center of mass 
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In generating the rigid partition using a symbolic manipulation software, consider 
the yoke body and the flexible domain as a single rigid body defined in the current 
configuration. 


Application of the virtual work principle to the flexible appendage yields 

; dV 


I 5^f (f-p£ F )dV = I 5e:o< 

Jvp Jv F 


(6.15) 


with 


£f = R + A + r F + T) s tJ R + F b A + F* ( rp + £ ) 


^F = Fb[ii + mbXu + iObxA + ©b x (fife x a) ] 

+ Fj [ ak x ( r F + n ) + 2 X £ + X (ojk X ( r F + £ )) + jj] = Fj a b = fJ ^ 


5^f = Fb [5x +6fib x a] + Fk [Sfibc x ( if + £ ) + 8n] 


where 


u = velocity of body frame wrt inertial frame 

iQb = angular rate of body frame wrt inertial frame 

coy = angular rate of appendage frame wrt inertial frame 

8x = virtual displacement of the body frame 

8fib = virtual rotation of the body frame 

5£k = virtual rotation of the appendage frame 

8n = virtual relative displacement 


The appendage angular rate can be decomposed into 

cok = F^ mk = w b + Wbk = Fb tOb + Fj iOb k ( 6 . 16 ) 

Similarly, 

50k = Ff 50k = FI 50b + Fj 80 bk (6.17) 

where 

Qbk = relative angular rate of appendage frame wrt body frame 
5Qbk = relative virtual rotation of the appendage frame wrt body frame 
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The frames are transformed with 


fb = Cbkfk 

Corresponding to Equation (6.9), the rigid degrees of freedom define 


(ur) = \ ffib j 
\ ffibk ) 

Substitution into Equation (6.15) yields following relations. 


Cbk H P dV 


W VF 


W VF 


[a x (Cbk ii) + Cbk (( u + n )x ii)] p dv 
If (if+ n)x ripdv) 


[B] {u F > 


/Vf 


CbkfrdV 




(a + Cbk ( if + n )) x (Cbk £f) dv 


r Vp 


( If + n ) x f F dv 


/vp 


2 Cbk X n) P dv 


r vp 


2 1 

J Vf 


2 AX (Cbk ( ffik x II ) ) + Cbk (( EF + n ) x ( fflk X n )) p dv 


(2 1 (( rp + n ) X ( X 11 )) p dV 


(6.18) 


(6.19) 


( 6 . 20 ) 


{Rrf} 


( 6 . 21 ) 
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The equations for the flexible domain results in 


I 

J Vf 


8n T f F dV 



2k P dV 


+ 



8e : adV 


( 6 . 22 ) 


Equations (6.20) - (6.22) yield necessary equations to generate the full equations of motion 
for a vehicle with articulated flexible appendage. All terms are retained. Vehicle specific 
truncation of nonlinear terms should be made to these equations. Detailed derivation using 
an alternate approach of a vehicle with articulated flexible appendage is provided in [5], 


Validation 


To validate the concept of mating flexible domain equations with code generated 
using a symbolic manipulation software, the HS-601 satellite [5] is chosen for 
implementation. AUTOLEV software generated the rigid body portion of the code. Rest 
of the code that deals with flexible domain has to be coded by the user. 

The specific example applies 40 ft-lbs of torque about the roll-axis for .1 sec. The 
magnitude of the torque characterizes the authority of the reaction control thrusters. The 
angular positions and rates are shown in Figures 6.3 - 6.4. The results are identical to the 
previous simulation that was implemented using the derivation described in [5]. Note that 
the response shown characterizes an actual vehicle. The effects of the flexibility are 
pronounced. In the rate plots, the magnitude of the rates peak at about 4 times the rate 
expected for a rigid vehicle with same applied torque. 

The use of a symbolic manipulation software with flexible domain equations offers 
large time saving in terms of both derivation and implementation. In this context, this 
approach is practical as an engineering tool. 
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7. Simulation Performance 

Digital simulations employing the hybrid dynamic model with the complete, 
nonlinear body -stabilized dynamics have been performed to demonstrate the stationkeeping 
control performance. Key simulation parameters which were added to produce worst case 
transient errors are: (1) 1.0 inches spacecraft CG offset from the pressure center of the 
south maneuver thrusters along the z-axis, (2) 5 Ibf thrusters with 5% thrust mismatch 
producing the worst case acceleration disturbance, (3) 1 .0 deg thruster misalignment in the 
direction adding disturbance, (4) flexible north and south wings at 0 deg wing angle, (5) 2 
RTI control transport delay, (6) solar wing drive with ±0.25 deg deadband, (7) thruster 
EPW (Electrical Impulse Width) error model, which computes the thruster impulse °n time 
delay as a function of time since last pulse. To accommodate the control loss due to E 
error, a fixed 4 msec thruster delay compensation was added to the command pulsewidth, 
(8) 8 msec pulsing constraint, (9) momentum wheel spinning at 45 ft-lb-sec throughout 
maneuver, (10) 0.035 deg sample to sample three sigma earth sensor noise, and (1 1) gyro 
sensor noises: rate random walk at PSD = (10‘ 7 d/s 2 ) 2 /Hz, angle random walk at PSD = 
(10- 4 d/s 2 ) 2 /Hz, angle noise at PSD = (1.3 x 10' 5 d/s 2 ) 2 /Hz and quantization of 0.3 
arcsec. 

With nominal structural mode frequencies and 0.5% structural damping, Figure 7.1 
shows the Stationkeeping Mode control performance from a 100-sec south maneuver. The 
spacecraft angular position and rate along with their estimates about each control axis are 
plotted. It also shows the acceleration bias estimate and the control acceleration command 
about each axis. The acceleration bias estimates in this run were initialized to zero The 
roll transient, which was induced primarily by the spacecraft yaw CG offset and the thrust 
mismatch, reaches 0.09 deg, while the yaw transient was primarily due to the thrust 
mismatch and was about 0.045 deg. The pitch transient induced by the combined effect ot 
canted and thrust mismatch is 0.01 deg. In the steady state, a limit cycle about the pitc 
axis resulted from the 8 msec thrusting constraint. The acceleration bias estimation 
converges within 5 sec, showing a smooth bias estimate in the steady state. Transient 
errors can be improved with an initialization of the acceleration bias estimates to their steady 
state values. Figure 7.2 shows the performance with both the roll and yaw acceleration 
bias estimates initialized to 4.5% off their steady state values. Due to an overestimate of 
the thrusting bias by 4.5%, the roll and yaw attitude were over controlled, yielding 
transients up to 0.06 deg and 0.042 deg, respectively, in the direction opposite to Figure 
7.1. Ideally, a perfect initialization would result in an significant reduction to the roll and 
yaw transients. 

Effects of structural mode uncertainty and damping were also investigated through 
simulations. Figure 7.3 shows the performance with a 100% frequency increase to each 
mode, and Figure 7.4 a 50% frequency drop. In both cases, a 0.25% structural damping 

was assumed. A 50% ( A ) frequency error is equivalent to a 75% [= 1-(1- A ) 2 ] change to 
stiffness of the flexible wings. It is not expected that the structural modes at the on orbit 
deployed condition would exceed more than 20% from the nominals. Likewise, the 
structural damping nominally exceeds 0.5%. The modal errors given above are to 
demonstrate the robustness of the control system. The simulation results show that 
variations in modal parameters have no major effects on the modal stability nor degrade the 
pointing performance. 
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100-SRC SOUTH MANEUVER 



Figure 7. 1 Stationkeeping Mode Control Performance From A 1 00 -sec 
South Maneuver with Zero Initial Acceleration Bias Estimates - Nominal 
Modal Frequencies, 0.5% Structural Damping 
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Figure 7.2 Stationkeeping Mode Control Performance From A 1 00 -sec Figure 7.1 (continued) 

South Maneuver with Acceleration Bias Estimates Initialized to 4.5% off 
Steady State Values - Nominal Modal Frequencies, 0.5% Structural Damping 

















8. Concluding Remarks 


The design and analysis of a stationkeeping control system for a body-stabilized 
spacecraft having flexible solar wings of 3 solar panels per wing were presented. The use 
of the hybrid coordinate modeling approach along with frequency domain analysis 
technique accurately modeled the rigid-flex coupling behavior. The design philosophy to 
stabilize the structural modes and to smooth the flexure was discussed The control system 
was designed to gain and/or phase-stabilize the structural modes. The lead inherent rate 
gyro references, the structural filters and the time-varying bias estimation gains were key 
factors to achieve a successful design. Control performance of Stationkeeping/Transition 
Modes during a south maneuver under the worst case simulation environment was 
demonstrated through digital simulation. The accuracy of the analytical model for structural 
mode/control loop interaction is best verified with the test data obtained from 
comprehensive ground testing. The control system as presented allows for a high degree 
of uncertainty on mode shape and frequency. 


An alternate, indirect implementation of flexible vehicle dynamics has been 
presented. The effort required to derive and to implement the equations of motion can be 
significantly reduced. A formulation through the virtual work principle allowed consistent 
derivation and discretization of the flexible domain within the context of the finite element 
method. Extensions into nonlinear flexible models can be made. 
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